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The chaotic hypothesis has several implications which have generated interest in the literature 
because of their generality and because a few exact predictions are among them. However its 
application to Physics problems requires attention and can lead to apparent inconsistencies. In 
particular there are several cases that have been considered in the literature in which singularities 
are built in the models: for instance when among the forces there are Lennard- Jones potentials 
(which are infinite in the origin) and the constraints imposed on the system do not forbid arbitrarily 
close approach to the singularity even though the average kinetic energy is bounded. The situation 
is well understood in certain special cases in which the system is subject to Gaussian noise; here 
the treatment of rather general singular systems is considered and the predictions of the chaotic 
hypothesis for such situations are derived. The main conclusion is that the chaotic hypothesis is 
perfectly adequate to describe the singular physical systems we consider, i.e. deterministic systems 
with thermostat forces acting according to Gauss' principle for the constraint of constant total kinetic 
energy ( "isokinetic Gaussian thermostats" ) , close and far from equilibrium. Near equilibrium it even 
predicts a fluctuation relation which, in deterministic cases with more general thermostat forces 
(i.e. not necessarily of Gaussian isokinetic nature), extends recent relations obtained in situations 
in which the thermostatting forces satisfy Gauss' principle. This relation agrees, where expected, 
with the fluctuation theorem for perfectly chaotic systems. The results are compared with some 
recent works in the literature. 
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Introduction 

There is a quite strong interest in stationary states of 
systems subject to the action of non conservative forces. 
These forces perform work on the system while by suit- 
able mechanisms heat is extracted, so that the system 
can stay in a statistically stationary state. Theoretical 
and experimental works are steadily becoming avalaible 
on the matter. Theoretical work implement the heat ex- 
traction in several ways introducing "thermostat mod- 
els" , which can be stochastic or deterministic forces. 

A strong idealization of a system in a nonequilibrium 
steady state subject to deterministic forces is provided 
by the Anosov systems: their motion can be considered 
to be paradigm of chaotic behavior, playing in chaotic 
dynamics the role that harmonic motions play in regu- 
lar dynamics. The chaoticity of the motions is imme- 
diately apparent from the definition of Anosov systems: 
locally around each point it has to be possible to draw 
three coordinate surfaces Ws, Wu, W such that segments 
of curves on Ws , Wu contract exponentially as time grows 
to +00 or, respectively, recedes to — oo while segments on 
W, the one dimensional "neutral" flow direction, neither 
expand nor contract. They change their length but keep 
it of the same order of the initial one. If the dynamics 
is described by a map the neutral direction is omitted in 
the definition. 

The chaotic hypothesis, [1, 2], see below, proposes that 
chaotic systems should be considered as Anosov systems 



"for practical purposes" . This has several consequences: 
in particular about fluctuations in time reversible models, 
where the hypothesis leads to severe constraints through 
the Fluctuation Theorem. This is a mathematical prop- 
erty of the large deviations function of the phase space 
contraction of a time reversible Anosov map S. How- 
ever, some obvious restrictions, analogous to the ones 
that are (often tacitly) assumed when one says that the 
"pendulum is isochronous" or that phonons in a crystal 
correspond to "harmonic excitations" , have to be taken 
into account when applying the hypothesis to realistic 
systems, which are not strictly Anosov systems. 

The prediction has been tested in several simulations 
and we summarize the precise statement of it below: 
usually the results have been positive. However, there 
have been, in the literature, a few claims of failure of the 
chaotic hypothesis based on the apparent failure of the 
predictions of the fluctuation theorem. Here we concen- 
trate on one such attempt, which studies systems vio- 
lating the Anosov property because singularities of the 
interparticle potentials play an important role in the dy- 
namics, [3] : a situation considered, correctly, in the liter- 
ature as not important for most physical properties but 
which requires care if the fluctuation relation is specifi- 
cally tested on such systems (in the same way care has 
to be used if isochrony is tested on a pendulum or har- 
monicity is tested in a crystal model). 

Here we show that even in singular systems the chaotic 
hypothesis and the fiuctuation theorem are not in con- 
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tradiction: wc develop a theory that extends the fluctu- 
ation theorem to singular systems continuing ideas that 
were introduced to study a special Gaussian noise ther- 
mostat, [4], (this is a "random thermostat" not to be 
confused with the thermostats satisfying Gauss' princi- 
ple for some non holonomic constraint, like the isokinetic 
constraint). 

The structure of the paper is the following: in section I 
we recall the basic notations and statements, and some 
alternative formulations of the fluctuation relation. We 
discuss its (trivial) form in equilibrium and how one can 
take a meaningful and non-trivial equilibrium limit. In 
section II we discuss the application of the chaotic hy- 
pothesis to singular systems. First we present a very 
simple example which shows that the effect of singulari- 
ties is very important. Then we discuss how one can ob- 
tain quantitative predictions on the modification of the 
fluctuation relation due to the presence of singularities. 
Finally we discuss a prescription to remove singularities 
that follows from a careful examination of the proof of 
the fluctuation theorem for Anosov flows. The results 
are compared with recent numerical simulations. In sec- 
tion III wc draw the conclusions and compare our inter- 
pretation with the one of [3]. 



I. THE FLUCTUATION RELATION 

We shall denote by 17 the phase space (a smooth com- 
pact boundaryless Riemannian manifold), by 5 : f2 — > f2 
an invertible map on fl and by a{x) the volume contrac- 
tion 



a{x) = — log I det dxS{x) 



(1) 



Time reversal is defined as an isometry I : fl< — >0 with 

IS = S-'^I, a{Ix) = -a{x) (2) 

If S is an Anosov maps, existence of a unique invariant 
probability distribution /i, called the SRB distribution 
and describing the long-time statistics of the motions 
whose initial data are chosen randomly with respect to 
the volume measure, is established, [5, 6]. It has the 
property that, with the exception of points x € f2 in a 
set of 0-volume, we have 



lim - Y,F{S*x)''^^{F)= [ FivUdy) 



(3) 



for all smooth observables F defined on phase space. 

It is intuitive that "phase space cannot expand"; this 
is expressed by the following result of Ruelle [7] : 

If (7+ (cr) it is cr+ > 

Clearly if 5 is volume preserving a+ = 0. If a+ > the 

system does not admit any stationary distribution of the 
form ii{dx) = p{x)dx, with density with respect to the 



volume mcasm-c dx (often called absolutely continuous 
with respect to the volume). 

This motivates calling systems for which (cr) > dis- 
sipative and conservative the others. 

For Anosov systems which are transitive {i.e. with a 
dense orbit), reversible and dissipative one can define the 
dimensionless phase space contraction, a quantity often 
related to entropy creation rate (see [8]), averaged over a 
time interval of size r. This is 



t/2-1 



(J_|_T ^ — ^ 

^ 7 _ /I 



(4) 



k=-T/2 



provided of course a+ > 0. 

Then for such systems the probability with respect to 

the stationary state, i.e. to the SRB distribution /i, 
that the variable p{x) takes values in A = \p,p + Sp] 
can be written as ni-(A) = g^^'^'^f^^ where 
C(p) is a suitable function and, for any fixed choice of 
A contained in an open interval {—p*,p*), p* > 1, the 
correction term at the exponent is 0(1) with respect to 
T~^, as T oo (this is often informally expressed as 
lim^^oo 7 logn^(p) = dp) for -p* < p < p*). The func- 
tion C(p) is called in probability theory the rate function 
for the large deviations of p. 

The function ({p) is analytic in p and convex in the in- 
terval of definition {—p*,p*). Analyticity and convexity 
of large deviation rates are general properties, established 
by Sinai and valid for the SRB-averages of smooth ob- 
servables (in Anosov systems), [6, 9, 10]. In fact more 
can be said for the specific case of the large deviation 
rate of the observable p, and one can prove the following 
fluctuation theorem: 

In transitive time reversible dissipative Anosov systems 
the rate function (^(p) for the dimensionless phase space 
contraction p{x) defined in (4) is analytic and strictly 
convex in an interval {—p*,p*) with +oo > p* > 1 
C{p) = "OO for \p\ > p* . Furthermore 



C{-p) = C{p) -p<y+, 



for 



\P\ < P 



(5) 



which is called the "fluctuation relation" (FR). 



Strict convexity follows from a theorem of Griffiths 
and Ruelle which shows that the only way strict con- 
vexity could fail is if 17(0;) = (p{Sx) — ip{x) + c where 
(p{x) is a smooth function (typically a Lipschitz continu- 
ous function) and c is a constant, see propositions (6.4.2) 
and (6.4.3) in [6]. The constant c vanishes if time re- 
versal holds and a-{x) = Lp(Sx) — (f{x) contradicts the 

assumption that cr+ > 0, because T^l^J^ (^{S''x) = 

[ifiiS^/'^-'^x) - (fiS-'/^x)] ^ as T ^ 00. 
The value of p* must be p* > 1 otherwise the average 
of p could not be 1 (as it is by its very definition) : it is 
defined, adopting the natural convention that ({p) = —00 
for the values of p whose probability goes to with t 
faster than exponentially, as the infimum of the p > for 
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which = —oo. Alternatively ±p* are the asymptotic 
slopes as A — > ±cxd of the Laplace transform log{e^P) srb 

The fluctuation relation was discovered in a numeri- 
cal experiment, [12], dealing with a non smooth system 
(hence not Anosov). The formulation and proof of the 
above proposition is in [1] and in the context of Anosov 
systems the relation (5) is properly called the fluctuation 
theorem. The difference between this theorem and other 
fluctuation relations proposed in the literature has been 
clarified in [13]. The theorem can be extended to Anosov 
fiows {i.e. to systems evolving in continuous time), [14]. 

Alternative formulations 

Sometimes, e.g. in [3, 15], rather than the above p 
the quantity a = J2^=~r/2 (^{S-'x) is considered and 
eq.(5) becomes 

C(— o) = C('^) ~ for I'^l < a* = p*cr+ (6) 

where (^{a) is trivially related to C{p). This form danger- 
ously suggests that in systems with (t+ = the distri- 
bution of a is asymmetric (because the extra condition 
\a\ < p*(J+ might be forgotten, see [16]). 

Note that p* is certainly < +oo because the variable 
a{x) is bounded (being continuous on the bounded man- 
ifold on which the Anosov map is defined). 

However no confusion should be made between p*(y+ 

and a,nax '==' max |cr(a;)|: unlike amax the quantity p* is a 
non trivial dynamical quantity, independent on the met- 
ric used on phase space to measure distances, hence vol- 
ume. This point has not been always understood and 
confusion has appeared in the published literature with 
unexpected consequences. In fact it is very easy to build 
examples of Anosov systems in which p*iT+ < (Jmax' still, 
this does not mean that fluctuation relation is violated 
for such systems. Some explicit examples are discussed 
in next Section. 



distribution. The distribution of a is symmetric (trivially 
by time reversal symmetry) and becomes a delta function 
around as r ^ cx). 

Nevertheless the fluctuation relation is non trivial in 
cases in which the map S depends on parameters E_ = 
{El, . . . ,En) and becomes volume preserving ("conserva- 
tive" ) as ^ ^ 0: in this case cr+ ^ as £ — > and one 
has to rewrite the fluctuation relation in an appropriate 
way to take a meaningful limit. 

The result is that the limit as £ —> of the fluctua- 
tion relation in which both sides are divided by makes 
sense and yields (in the case considered here of transitive 
Anosov dynamical systems) relations which are non triv- 
ial and that can be interpreted as giving Green-Kubo for- 
mulae and Onsager reciprocity for transport coeflacients, 
[17, 18]. 

In fact the very definition of the duality between cur- 
rents and fluxes so familiar in nonequilibrium thermo- 
dynamics since Onsager can be set up in such systems 
using as generating function the a+ regarded as a func- 
tion of E_. Note that the fluxes are usually "currents" 
divided by the temperature: therefore via the above in- 
terpretation one can try to deflne the temperature even 
in nonequilibrium situations, [8, 19, 20]. 



II. SINGULAR SYSTEMS 

The fluctuation relation has been proved only for 

Anosov systems. However, a Chaotic Hypothesis has been 
proposed, which states that, for the purpose of study- 
ing the physically interesting observables, a chaotic dy- 
namical system can be considered as an Anosov system, 
[1, 21, 22]. 

In applying the chaotic hypothesis to singular systems, 
e.g. a system of particles interacting via a Lennard- Jones 
potential (which is infinite in the origin), one might en- 
counter apparent difficulties. We will discuss them in the 
following. 



Conservative systems and the equilibrium limit 



The effect of a change of metric for Anosov flows 



Considering more closely the cases a+ = it follows 
that a{x) = ip{Sx) — if{x) (again by the above mentioned 
result of Griffiths and Ruelle), with </? a smooth function 
of phase space. Hence the variable 



r/2-l 



a = - ^ g{S^x) = 



_ <f{S ^x) — ip{S^x) 



(7) 



j=-T/2 



is bounded and tends to uniformly. One could repeat 
the theory developed for p when c7_)_ > but one would 
reach the conclusion that ({a) = — oo for |a| > and we 
sec that the result is trivial. In fact in this case it follows 
that the system admits an absolutely continuous SRB 



The simplest example (out of many) is provided by the 
simplest conservative system which is strictly an Anosov 
transitive system and which has therefore an SRB dis- 
tribution: this is the geodesic flow St on a surface of 
constant negative curvature, [23]. We discuss here an 
evolution in continuous time because the matter is con- 
sidered in the literature for such systems, [16] (even sim- 
pler examples are possible for time evolution maps). 

The phase space M is compact, time reversal is just 
momentum reversal and the natural metric, induced by 
the Lobatchevsky metric gij{q) on the surface, is time 
reversal invariant: the SRB distribution is the Liouville 
distribution and ^{x) = 0. However one can introduce 
a function ^{x) on M which is very large in a small 
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vicinity of a point xq, arbitrarily selected, constant out- 
side a slightly larger vicinity of xq and positive every- 
where. A new metric could be defined as gnew{x) = 

{^{x) + ^{Ix))g{x) e^^^'^hj{x): it is still time rever- 
sal invariant but its volume elements will no longer be 
invariant under the time evolution St associated with 
the geodesic flow with respect to the Lobatchcvsky met- 
ric. The rate of change of phase space volume in the new 
metric will be o'new 

[x) = ^F{x). Then the phase space 

contraction (Jnew{x) takes values that not only are not 
identically but which can in general be arbitrarily large, 
depending on the specific choice of The distribu- 

tion of a = i (Jnew{StX)dt = T~'^[F{SrX) - F{x)], 

at any finite time, will violate (6), simply because it is 
symmetric around 0, by time reversal. 

In the limit r — > oo, as long as F(x) is bounded, a = 
T-^[F{Srx) - F{x)] -f^55*0 uniformly in X, as in the 
corresponding map case, and the SRB distribution of a 
will tend to a delta function centered in (hence C(a) = 
—00 for a 7^ 0). However, if F{x) is not bounded {e.g. if 
it is allowed to become infinite in xq) this is not the case 
in general, as we shall discuss in detail in next section. 



The effect of singular boundeiry terms 

One can realize that terms of the form [F{Srx) — 
F(x)] with F{x) not bounded can affect the large fluctu- 
ations of (t{x), at least if the probability of an arbitrarily 
large value of F is not too small, i.e. if asymptotically 
for big values of F it is exponentially small in F (or 
larger), e.g. it is of the form ~ e"**^, for some constant 
K > 0. This is a valuable and interesting remark brought 
up for the first time, and correctly interpreted, already 
in [4] and in the following papers [24, 25]. The analysis 
of [4, 24, 25] applies to cases where the unbounded fluc- 
tuations are driven by an external white noise. In the 
following we extend the theoretical analysis in [24, 25] to 
cases in which the unbounded fluctuations do not arise 
from a Gaussian noise but from a deterministic evolution 
like the ones in [3, 15]: this is a simple extension of the 
main idea and method of [24] and provides an alternative 
interpretation to the analysis in [3, 15]. 

Our analysis can be applied to the example of the 
Anosov flow with singular metric considered above and to 
more realistic systems: among them systems of particles 
interacting via an unbounded potential (like a Lennard- 
Jones (LJ) or a Weeks-Chandler-Andersen (WCA) po- 
tential), driven by an external field and subject to an 
isokinetic or a Nose Hoover thermostat. To be definite 
one can consider system of A'^ particles in d dimensions, 
described by evolution equations Pj = E — 9q;$ — api, 
cji = Pi- For an isokinetic Gaussian thermostat, a is a 
function of Pi, chosen so to keep the total kinetic energy 
fixed to p| = Nd(3~^. For a Nose-Hoover thermostat 
a{t) is a variable independent of (li{t),Pi{t) and satisfy 



Q, /? > parameters. 

In both cases the phase space contraction a{x) has the 
form ao{x) — (3^V{x), where f3 has the interpretation 
of inverse temperature. In the isokinetic case, a^o^x) is 
bounded, and y = $. In the Nose-Hoover case (To{x) 
has, in the SRB distribution, a fast decaying tail (Gaus- 
sian at equilibrium, and likely to remain such in presence 

of external forcing) and V 4+*(q)+'5^' [26, 27]. 

In both cases, in equilibrium, the SRB probability of V 
has an exponential tail ~ e~^'^ (possibly with power-law 
corrections). For the purpose of illustration we assume, 
from now on, that the same happens in presence of the 
force E. This is an essential and far from obvious as- 
sumption useful, as discussed below, to understand the 
possible role of the singularities, but it should not be 
assumed lightly as it is well known that the SRB distri- 
butions may have very peculiar E dependence and, at 
the moment, a not intuitive character, [28, 29]. Nev- 
ertheless, in preliminary numerical simulations, it seems 
approximately correct, at least within the accuracy of the 
numerical data and for [E[ not too large; furthermore the 
analysis that follows can be naturally adapted to more 
general assumptions on the tails. 

In such cases the non normalized variable a, (introduced 
before Eq. (6)) has the form ao+^iVi — Vf) where Vi, Vf 
are the values of V{x) at the initial and final instants 
of the time interval of size r on which a is defined, and 

ao'^= ^ Jq aoiStx)dt: 



a=- [ a{Stx)dt = ao + - (Vi - Vf) 



(8) 



If the system is chaotic and r is large, the variables 
<io,^)^/ can be regarded as independently distributed, 
because ao depends essentially only on the length r of 
the time interval, while Vi and Vf depend on the precise 
locations of the extremes of the interval. Moreover the 
distribution ofV = ViOvV = Vf ]s essentially ~ e~'^^dV 
to leading order as V ^ oc, as discussed above. There- 
fore the rate function of the variable a can be computed 
as 



lim — log / da, 

. pr(o(ao)-0Vi-0Vf> 



/ dV / 
Jo Jo 



dVf 



'S[T{a-ao) + pVi-pVf] 
= lim - log r ^ daoe^^''(''°'-^l''-"''l 

T^oo T J-p'a+ 

where Co(ao) is the rate function of ao; thus 



(9) 



C(a) 



max 

aoe[-p*(T+,p*<7+] 



Co(ao) - [a - ao\ 



(10) 



±1, by the strict convexity of 



ing the evolution equation a = i [ J2i Pi — Nd(3 ^] , with 



Defining by CoIot) 
Co(ao) it follows 

Co(a_) — a_+a , a < a_ 

C(a) = { Co(a) , ae [a_,a+] (11) 

Co(a-i-) + a-i- - a , a > a+ 
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If we assume that Co(ao) satisfies FR (as expected 
from the chaotic hypothesis, see below), then Co(ao) — 
Co(~cio) + o-o SLud by differentiation it follows that a_ = 
— 0"+, where (t+ is the location of the maximum of Coj 
i.e. is the average of a, and that Co(q.-) = Co(^o'+) = 
Co(o'+) ~ cr+ = ~cr+. Moreover it is clear that a+ > (t+ 
because Co(fl+) < 0. Using these informations one can 
show that, for a > 0: 

{a , a < (7+ 

Co(a) + a , (T+<a<a+ (12) 
Co(a+) + a+ , a > a+ 

It follows that, if Co(ao) satisfies FR up to a — p*aj^, then 
C(a) satisfies FR only in the interval \a\ < |a_| = cr+. 
Outside this interval ({a) does not satisfy the FR and 
in particular for a > a+ it is C(a) — C(~a) = const., as 
already described in [24]. Eq. (12) is the generalization of 
the result of [24] to the case where Co(ao) is not Gaussian. 

Translated into the normalized variables po = aQ/(T+ 
and p = a/a+, this means that, even if the rate func- 
tion of Po satisfies FR up to p* > 1, the rate function 
of p verifies FR only for \p\ < 1. This is the effect due 
to the presence of the singular boundary term. Note 
that the scenario above applies only to the case in which 
Vi, Vf are unbounded and have exponential tails. A rep- 
etition of the discussion above in the case that Vi, Vf are 
unbounded but with tails faster than exponential would 
lead to the conclusion that C,[a) ~ Co (a)- In particular 
if Vi,Vf are assumed to be bounded ({a) — Co (a)- Of 
course in these cases the times of convergence of C(a) 
to Co (a) will depend on the details of the tails of Vi,Vf 
(for instance if Vi, Vf are bounded by a constant B, the 
times of convergence will grow with B). Note also that 
the result above does not depend on the details of the 
distribution of Vi,Vf for small V (in particular it does 
not depend on the lower cutoff V = assumed in Eq. 9). 

An example of C(a) is reported in Fig. 1: it is a sim- 
ple stochastic model for the FT (taken from Sect. 5 
in [30], see also the extensions in [31, 32]). The ex- 
ample is the Ising model without interaction in a field 
h, i.e. a Bernoulli scheme with symbols ± with prob- 
abilities p± = 2^5or- Defining oq = 7 LUo^ ^/iCTi, so 

that CT+ = (oo) = 2/itanh/i, and setting x'^= l±2i>ll2hl^ 
and s{x) = — xlogx — (1 — x) log(l — x), one computes 
Co(ao) = s{x) + iflo + const which is not Gaussian and it 
is defined in the interval [—a*, a*] with a* — 2h. In this 
case the large deviation function Co(ao) satisfies FR for 

|ao| <a*. If a singular term = -log(X;j=o 2"'"^^) 
is added to oq, defining a = ao + l3{Vi — Vf) (with 
f3 = log2(l -t- e^'*) so that the probability distribution 
of 1^ is e~^^ for large V), the resulting C(a) does not 
verify FR for a > (a) — 2/itanh/i. In particular, for 
h ^ 0, the interval in which the FR is satisfied vanishes. 




FIG. 1: An example in a stochastic model of FR. The graph 
gives the two functions Co(i) and ({a) for h = 0., 0.25, 0.5. 
The average of a is (a) = (j+ = 2/itanh/i, a+ = 2fttanh3/i 
and a* = 2h. The function ("(a) is obtained from Co(fl) by 
continuing it for a < a_ = — (t+ and a > a+ with straight 
lines of slope ±1. It does not satisfy the FR for jaj > (a). 
As /i ^ 0, (a) 0, which means that the interval in which 
the FR is verified shrinks to 0. In this limit a+ — > 0, so C^(a) 
approaches —\a\ (dashed lines). Rephrasing this in terms of 
p — one obtains that FR remains always valid for \p\ < 1, 

even as /i — > 0. The three curves for Co(ffl) have the same 
tangent on left side. The function Co (a) is finite only in the 
interval [—2h,2h] and it is —00 outside it, while the function 
("(a) is finite for all a'a and is a straight line outside [a_,a+]. 



How to remove singularities 

From the discussion above it turns out that singular 
terms which are proportional to total derivatives of un- 
bounded functions (like the term ^ that appears in the 
phase space contraction rate of thermostatted systems) 
can induce "undesired" (or "unphysical" ) modifications 
of the large deviations function C(p)- 

On heuristic grounds, when dealing with singular sys- 
tems, one could follow the prescription that unbounded 
terms in a{x) which are proportional to total derivatives 
should be subtracted from the phase space contraction 
rate. If the resulting 170(2;) is bounded (as it is e.g. for 
the Gaussian isokinetic thermostat models considered) or 
at least if the tails of its distribution decay faster than 
exponentially, then its large deviations function should 
verify the FR for \p\ < p* , p* being the intrinsic dynamic 
quantity defined above. 

Note that after the subtraction of the divergent terms 
the remaining contraction, in the considered cases, is 
bounded for isokinetic thermostats or has a tail decay- 
ing faster than exponential in the case of Nose-Hoover 
thermostats. In the following for definiteness we will as- 
sume ao bounded but the same discussion is valid for ao 
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unbounded with tails decaying faster than exponcntiaL 

If the singular terms are not subtracted, the FR will 
appear to be valid only for \p\ < 1 even if p* > 1. This 
seems to have generated statements that the Chaotic Hy- 
pothesis does not apply to isokinetic systems, see [3]. 

The heuristic prescription above can be motivated by 
a careful analysis of the proof of the fluctuation theo- 
rem for Anosov flows. In the following let us call again 
a the integral of the total phase space contraction rate 
a{x) (which includes singular terms) and ao the integral 
of the bounded variable aQ{x) from which singular total 
derivatives have been removed. 

The fluctuation theorem was proved in [1, 2, 33] for 
Anosov maps and only later it has been extended in [14] 
to Anosov flows. Very sketchily, the extension of the 
fluctuation theorem to Anosov flows in [14] is proved as 
follows. One reduces the Anosov flow on f2 to a map 
via a Poincarc's section, associated with surfaces on fl 
transversal to the flow. The passage of the flow through 
any one of such surfaces is called a timing event. The 
map between two consecutive timing events is called a 
"Poincarc's map" . The union Qp of the surfaces repre- 
sents the phase space of the Poincarc's map. The sur- 
faces in fip can be suitably chosen, in such a way that 
the Poincare's map is a chaotic map which although not 
smooth, hence not an Anosov map, has (a non trivial fact 
[14]), all the properties necessary to prove the fluctuation 
theorem (which therefore applies to systems more general 
than the Anosov maps, although there is not a general 
characterization of the systems which are not Anosov and 
to which it applies). So, for such a map the fluctuation 
theorem holds and this in turn leads to a FR for the 
flow by the theory in [14] under the assumption that the 
variable a{x) is bounded. 

If, as in the case under analysis, (t(.t) is not bounded, 
we can interpret the chaotic hypothesis as applying to the 
map associated with a Poincare's section which avoids 
the singularities of the potential, a very natural prescrip- 
tion which allows us to apply the theory in [14] and derive 
a FR for both the map and the flow. For instance, we can 
choose as timing events the instants in which either the 
potential energy or the Nose's "extended Hamiltonian" 
exceed some fixed value V. If we make this choice, the 
(discrete) average a of a{x) over a sequence of iterations 
of the Poincare's map will coincide with the (discrete) 
average oq of <7q(x) along the same sequence: this simply 
follows from the remark that by construction the total in- 
crement of (t{x) — (To (x) between two timing events, given 
by l3{Vf — Vi), is (by construction Qp is chosen as a sub- 
set of {x e n : V{x) = V} where Vf = Vi). Then, by the 
same argument in [14], the fact that the rate function of 
ao satisfies a FR and that cro(x) is bounded implies that 
the rate function of the continuous average ao of ao{x) 
along a trajectory of the flow will satisfy the fluctuation 
theorem. 

Therefore the distribution of ao will satisfy the FR {by 

the chaotic hypothesis) for ]ao| < p*o'-^-. By the above 
maximum argument, the distribution of a will also verify, 



as a consequence, the FR but only for ]a| < cr+, i.e. in 
the form (12). 

Then the (natural) prescription to study FR for chaotic 
flows is to reduce the problem to a chaotic map consid- 
ering only Poincare's sections which do not pass through 
a singularity of a{x). The sum of cr{x) over a large num- 
ber of timing events on such sections is equal to the time 
integral of ao{x) plus a bounded term which can be ne- 
glected. Thus the prescription on the choice of Poincare's 
sections is equivalent to the heuristic prescription of re- 
moving from cr(x) all the unbounded total derivatives. 

It follows that the chaotic hypothesis leads to a clear 
prediction on the outcome of possible numerical simula- 
tions of particle systems interacting via unbounded po- 
tentials and subject to the isokinetic or the Nose-Hoover 
thermostat: the FR will hold for all ]a| < cr+ and, once 
the term ^ is removed, for all [ao] < P*<t+ with p* > 1. 
Note that in the cases under analysis ao coincides with 
the dissipation function of Evans and Scarles that was in 
fact predicted to satisfy FR [3, 15], even though for dif- 
ferent reasons. We believe that the correct interpretation 
of the fact that FR for Coi'^o) holds for all [ao] < p*(J+ is 
the one given above. 

The numerical results of [34-36] agree with the pre- 
diction that FR for the rate function of ao is valid even 
beyond ag = (7+ . The prediction that (at least near equi- 
librium) the rate function of a should satisfy FR only up 
to a = cr+ and that should become linear for a > a+ 
at the moment has been experimentally confirmed only 
in Gaussian cases [4, 24, 25]. It would be very interest- 
ing to investigate in detail the structure of C(a) even in 
non Gaussian cases. Note that this is far from being an 
easy task (in particular the analysis in [36] was not so- 
phisticated enough to study this problem). In fact, as 
discussed in detail in [35], the presence in the definition 
of (T of a total derivative of an unbounded function may 
enlarge of 2 orders of magnitudes the times needed for 
the probability distribution of a to reach its asymptotic 
shape: even in the Gaussian region (small fluctuations of 
a around cr+) the convergence times for C(a) are found 
to be of order 1000 decorrelation times, versus a time 
of order 10 decorrelation times needed for Co(flo) to con- 
verge to its asymptotic shape [35]. Clearly, for times of 
order 1000 decorrelation times, it is very hard to observe 
fluctuations of a larger than a+ — (t+ . In order to verify 
the prediction for the shape of ({a) beyond a = a+ an 
experiment speciflcally designed for this purpose would 
be needed, together with a detailed investigation of the 
flnite time corrections to C{a), along the lines in [36]. 



III. CONCLUSIONS AND REMARKS 

We showed that the Chaotic Hypothesis can be ap- 
plied even to singular chaotic systems (in particular even 

to Gaussian isokinetic or Nose Hoover thermostatted 
systems), by identifying their macroscopic behavior 
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with that of reversible Anosov systems with singular 
metric. Reversible Anosov systems with singular metric 
are systems to which the mathematical analysis usually 
leading to FR can be (rigorously) repeated to lead to a 
modified FR, illustrated by Eq.s (11), (12). Note that for 
(7+ — > Eq. (11) tends to the distribution ({a) = —\a\ 
violating the usual FR (simply because the limiting 
distribution is symmetric, by time reversal). For Anosov 
systems with singular metric, the prescription to avoid 
oddities {i.e. to avoid a modified FR) is to subtract from 
a{x) a total derivative dV/dt, in such a way that the 
variable (Tq = c — dV/dt is bounded or has faster than 
exponential tails. The distribution of ctq will verify the 
FR also for \p\ > 1. This prescription is equivalent to the 
very reasonable prescription that the Poincarc's section 
used for mapping the fiow into a map does not pass 
through a singularity of (t{x). Accepting the Chaotic 
Hypothesis, we propose to apply the same prescription 
to remove singularities to singular chaotic systems. Our 
prescription coincides with other prescriptions proposed 
earlier (for different reasons) in the literature. 
The analysis in Sect. 2 above applies as well to under- 
stand how to apply the FR to systems with Gaussian (or 
unbounded) noise and the compatibility between the gen- 
eral theory of [31, 32, 37] with the works [24] and [25, 38]. 

Note that the picture we propose is different from 
the interpretation of the apparent violations to FR 
in singular systems proposed recently in [3], where in 
particular it is argued that FR and CH do not hold for 
thermostatted systems near equilibrium. We conclude 
by comparing more closely our discussion with the 
corresponding discussion in [3]. 

(1) As stressed above it is dangerous (and wrong) 
to consider (6) without the restriction \a\ < p*a+ as the 
prediction of fluctuation theorem. In [3] the authors, 
after having correctly pointed out this point, seem (quite 
surprisingly!) to forget about this condition in the 
following. For instance, when studying the problem of 
approach to equilibrium, in order to show a contradiction 
between FR and GK relations, they assume that the 
relation Eq.(6) without the condition \a\ < p*a+ "is 
correct both at equilibrium and near equilibrium" and 
they proceed to infer from this a contradiction. Of 
course such assumption is wrong and the fact that from 
this contradictions follow is not an argument against CH 
or FR. 

(2) An argument in [3] is supposed to prove that the rela- 
tion C{—a) = ({a) — a (without the condition \a\ < P*<t+) 
holds for reversible Anosov systems for all a's, also for 
cr+ = (in particular they say that "i/ie division by (T+ 
does not seem to be necessary for the proof in /57/'). 
This is not the case: at equilibrium as well as near equi- 



librium, as remarked above and as illustrated also by 
[23], there are examples of systems for which the proof 
of FR can be rigorously repeated step by step but for 
which the correct^ conclusion of the proof is that the re- 
lation Ci~C') = C('^) — a is violated for a > p*(J+. For 
instance, this is the case for a conservative Anosov flow 
with singular metric (in which the relation above is vio- 
lated trivially by time reversal). These counterexamples 
show that the assumption (T+ > is. instead, essential 
for the proof of fluctuation theorem. The necessity of 
the assumption fT+ > is stressed in the early paper [11] 
which the Authors of [3] quote; it is stressed also in the 

paper [2] which also makes clear that C(~^) = C(<^) ~ o, 
can only hold under the assumption that |a| does not 

exceed a maximum vahic. 

(3) The analysis in Sect. 2 above shows that the prob- 
ability distribution describing isokinetic systems near 
equilibrium are SRB distributions (contrary to what is 
claimed in [3]): this is mathematically obvious by the 
very definition of SRB distribution in the case of Anosov 
systems {even if isokinetic or in general with singular 
metric, see e.g. the geodesic flow discussed above) and it 
appears to be true also in non Anosov systems that have 
so far been considered. 

(4) In the case of the thermostatted particle systems con- 
sidered in [3] the unbounded derivative ^ is also the 
contraction rate of the volume in equilibrium, i.e. with 
E = 0. Thus, for E 0, one can remove the total deriva- 
tive from (t(.t) simply considering the contraction with 
respect to the equilibrium invariant distribution e~^^ , as 
stated in [3, 15]. However, this observation docs not pro- 
vide a general prescription to remove the singular part 
from the phase space contraction rate because it rests on 
the very special fact that the singularities of the function 
V{x) {i.e. of the potential $) do not depend on E. The 
prescription that the phase space contraction should be 
computed on non singular Poincarc's sections, instead, 
does not require any other assumption. In general the 
two prescriptions and the corresponding predictions differ 
and we believe that in general the prescription of comput- 
ing a with respect to the eqiiilibrium invariant distribu- 
tion has not the desired effect of removing all singularities 
(then in general cr with respect to the equilibrium invari- 
ant distribution could violate FR for 0-+ < a < p*a+). 
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